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Abstract 

We explain the linear growth of smooth solid helium facets by the presence of lattice point defects. To 
implement this task, the framework of very general two- velocity elasticity theory equations is developed. 
Boundary conditions for these equations for various surface types are derived. We also suggest additional 
experiments to justify the concept. 

PACS: 67.40.Bz, 67.40.Pm, 68.35. Ct, 61.72.Ji 

1 Introduction 

Existence of two distinct states of a crystal surface is well known: it may be either smooth or rough (for a 
review, see [0). A smooth surface is characterized by a long-range order and small fluctuations. On the 
contrary, a rough surface behaves differently — it does not exhibit a long-range order and its displacement 
fluctuates heavily. These equilibrium properties lead to different kinetic properties. While the rough surface 
is usually supposed to grow easily (as described by the growth coefficient), the smooth one (providing the 
crystal has no dislocations) is characterized by zero growth coefficient and grows with nuclei of new atomic 
layer. In accordance with this mechanism, one should not observe the linear growth rate if low overpressure 
is applied. Reality is different — experiments 0] demonstrate that a smooth helium surface free of screw 
dislocations grows linearly. This work is an attempt to explain this behavior by the presence of lattice point 
defects (vacancies). The idea is similar to that suggested by Herring || and by Lifshitz Q as an explanation 
of polycrystal flow. It is quite simple: the mass flux in bulk helium is attributed to the motion of vacancies. 
This flux is the mass transfer through the lattice. Therefore, if vacancies are allowed to be created on the 
bottom edge of the sample (the boundary between the crystal and the wall, see Fig. |l|) and to annihilate on 
the top of it (on the smooth crystal- liquid interface), then the crystal grows. 




silica wall 



Figure 1: Typical experimental layout 

The suggested crystal growth mechanism may be explained as follows. Since the smooth crystal facet (the 
top one in Fig. |l|) cannot move with respect to the lattice, it moves upward stuck to the lattice. Vacuities 
to appear due to this on the bottom edge of the crystal transform into lattice defects (ordinary vacancies) 
and go up through (and faster than) the bulk helium. They finally vanish in the liquid on the top smooth 
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surface of the sample. In other words, the crystal grows on the boundary between helium and the wall, 
rather than on the smooth solid-liquid interface (which nevertheless provides mass supply for the growth). 
It is important to emphasize that this scenario may occur if and only if the vacancies are allowed to emerge 
on the bottom edge of the crystal. One may say that this boundary is in some sense "atomically rough" - 
it can grow new atomic layers. For this condition to be satisfied, the wall surface should have a disordered 
shape or be slightly tilted with respect to basal planes of the crystal (thereby forming a vicinal interface). 
This ensures that the surface can play the role of a source or a sink of vacancies. An atomically flat wall 
parallel to the basal plane should, in contrast, behave like a normal smooth surface — it is fixed to the 
lattice. This is due to the fact that for the surface to move new atomic layer nuclei have to be created. 

This paper is organized as follows. In Sees. |^ and ^| we derive very general two- velocity elasticity theory 
equations. They consist of conditions for the conventional elasticity theory variables (including lattice 
velocity) and equations for a macroscopic description of the quasiparticle gas (including the quasiparticle 
gas velocity). 

Equations to be derived are similar to those of the two- velocity superfluid hydrodynamics. Velocities of 
lattice and excitation gas in our equations replace superfluid and normal component velocities of the two- 
fluid hydrodynamics. Similarly to the usual linear phonon hydrodynamics (see [|| § 71) Umklapp-process 
(which result in the nonconservation of the total quasimomentum in quasiparticle collisions) probability is 
supposed to be low. In the low-temperature region considered here this assumption is quite reasonable. 
We also neglect dissipation here. This means that our analysis is limited to the terms of the first order in 
gradients. 

Derivation procedure of exact (nonlinear) hydrodynamics equations for superfluid can be implemented 
(see H) from phenomenological considerations, using conservation laws. A constitutive argument for this 
derivation is the statement that the superfluid flow is potential. This is an intrinsic property of the order 
parameter in a superfluid. Such condition is unavailable for a crystal (moreover, there is no quasimomentum 
conservation relation in the nonlinear description, see Eq. ( |l5|) below). 

We deduce the two-velocity elasticity theory equations using a more general approach (see the paper by 
Pushkarov and one of the authors 0, as well as || and ||). It is based on the kinetic equation description of 
the quasiparticle dynamics. The realization of this technique, particularly in nonlinear situation, is a matter 
of considerable interest for its own sake not only for a solid but also for a superfluid (that this procedure 
is possible is mentioned in Q §77). With this technique, we find exact expressions for all hydrodynamic 
variables and their dependence (in terms of the quasiparticle energy spectrum) on the relative velocity of 
components. It is trivial to extend the equations obtained for the solid dynamics to the simpler case of 
superfluid hydrodynamics. 

Boundary conditions for our equations depend on the surface type; in Sec.^J we thoroughly consider three 
possibilities: rough (Sec. |4.1.l| ) and smooth (Sec. 4.1.2) interfaces between solid and liquid helium, and the 



rough boundary between solid helium and normal hard wall (Sec. 4.2.1). Finally, in Sec. || we calculate the 
growth rate for the crystal. 



2 Definitions 

Following the principles in Refs.0, ||, ^| we employ the Euler approach to the lattice description. We thus 
introduce three "node numbers" N a (a = 1,2,3). They are functions of space coordinates r and time t, 
N a = N a (r, t). From now on, Greek indices (like a here) are used for the "lattice space" and Latin indices 
{e.g., i in Xi for the components of r) for the real space components. Defining the reciprocal lattice vectors 
as a a = dN a /dr, we get the elementary lattice translation vectors as a Q a^ = 6^. Taking the time 

derivative, we obtain the lattice velocity as w = — a a N a . The elastic energy Ei of the lattice is a function 
of the deformation. Moreover, since it depends not on the spatial orientation of the infinitesimal sample 
(the space is isotropic), but on the relative position of the af vectors, we may write Ei — Ei (g Q/3 ), where 
g a @ = a Q a /3 is a symmetric "metric tensor" of the lattice space. 

We are now ready to describe quasiparticle degrees of freedom. We do not specify the quasiparticle 
nature at the moment (be it phonons, vacancies as in [^J, or electrons as in Q). Actually, all the equations 
below imply the summation over all branches of excitations; we do not explicitly write the sum for brevity. 
Any quasiparticle should be characterized by its mass m (zero for phonons, positive for electrons, and 
negative for vacancies), coordinate, and momentum. Since quasiparticles exist in the lattice background, the 
quasimomentum should be used. Their energy in the frame of reference of the lattice e = e(a Q (p — mw), g a °) 
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is a periodic function of the quasimomentum p (its periods are 2nUa^). In laboratory frame of reference, we 
have the quasiparticle energy (see Q]) 

de w 2 

e = e + tow— h to — . 

op 2 

We also use the variables k = p — tow and k a — a Q k. Quasiparticle dynamics is determined by the Hamilton 
function 

H = e + pw — mw 2 /2. 

We now introduce the distribution function /(r, p) (it is also a periodic function of the quasimomentum p). 
Its kinetics is governed by the Boltzmann equation 

df_ + d£dH_ _ df_dH_ _ 
dt dr dp dp dr 

Using this distribution function, we can obtain macroscopic quantities like the mass density as 

1/2 

p = pi + (mf) = M dot (g al3 ) + run, 

where the angle brackets denote the integration over quasimomentum space, {) = J d 3 p/(27r?i) 3 , pi is the 
lattice density, M is the mass of an elementary cell, and n = (/). 

We consider a quasi-equilibrium distribution function. The complete set of quantities conserved in quasi- 
particle collisions consists of their mass (proportional to their quantity for "real" particles like electrons 
and vacancies and zero for phonons), energy, and quasimomentum (in the low-temperature region, Umklapp 
processes may be neglected). Consequently, the most general quasi-equilibrium distribution is a function of 

e — kv — m/io e — (p — mw)v — mpo e — pv — to/^q + towv e — pv — m<j) 



rj~i rj~i rj^i rj~i 

where (iq — wv = <f>. The Lagrange coefficients T, v, and po denote temperature, the velocity relative to the 
lattice, and chemical potential of the quasiparticle gas, accordingly. For definiteness, we assume that the 
excitations are Bose particles. The distribution function is then given by 

and ln((/ + 1)//) = z.Q 

We can now calculate other macroscopic parameters with this distribution function. For the mass flux, 
we have 

J = pi w + wtott- + j = pw + j = pw + TO71V, (3) 

where the mass flux with reference to the lattice is 

/.de\ 
j = to (/-) = mnv. 

Using J, we can write the mass conservation as 

P + Al = 0. 

The number of real (massive) particles is also conserved in the bulk, and we therefore have one additional 
conservation law 

pi + {piWi) t i = 0. 



x For further convenience, we also provide here the result of the distribution function integration: 

e z - 1 

/ dz = In = - ln(/ + 1). 
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Similarly, the energy density is given by 

E^ Pl ^ + E l (g"P) + (if) =P^ + E t (g°0) + wj + (ef) . (4) 

This equation allows us to prove (and find) exact macroscopic equivalents of the microscopic quantities 
introduced above. The total energy density of the crystal can be obtained via a Galilean transformation, 

E = E + ^f+3w, (5) 

where Eo — Eo(a a , S, p, K) is the energy in the frame of reference of the lattice, with K = (k/) characterizing 
the quasimomentum density. A reasonable expression for the Eq differential 

(IEq = \ija U j da" + T dS + fi dp + v dK (6) 

we obtain with the conventional definition of the entropy density for the Bose gas, 

S = ((/ + 1) ln(/ + 1) - / In/) = (fx + ln(/ + 1)) . 

Its differential being 

dS = ((ln(/ + 1) + 1 - In / - l)d/) = (zdf) . 
Subtracting the differentials of (Q) and (||), we obtain 

= A, ,</,., d< + TdS + pdp + vdK - dEi (g al} ) - d (ef) 

= \ ijauj d< + T (z df) + fi dp + v dK - dE t (g al3 ) - (e df) - (/ de) 

= X ijavj d< + ((Tz - e + vk) df) + //dp - dEi (g al3 ) + v (/ dk) - (/ de) (7) 
We now transform the part of this equation related to the lattice deformation 

%pig a p dg a ? - dEi + v (/ dk) - (/ de) = ^p l9af3 dg a ^ - ^ dg a ? - mv (/ dw) 

- { f [wk ((p - mw) da « - wa « dw) + ka d ^) ) 

dEi I ( ( de \ \\ A a0 , ( x / f d(Tx + pv) 



/ dEi J ( de \ \\ . 



+ / -53^ 2af<a^d< 



+ (Sij(T (ln(/ + 1)) + fipi) + ViPj - mnViWj)a„j da" 

w 2 \ 
TS + p— + Ei — E + Pv + ppi + mnpo J — Ay + ViPj — mnviWj I a„j da". 



Finally, from (]?]) we get: 
= (jj, — no)m dn 



+ | Ay - Ay + Sij I TS + Ei - E + pp + Pv + p— J + UjPj — mnv l Wj j a„j da" , (8) 
where we introduced Pi = (pif) and 

The terms in (^) are independent, and each of them must therefore be equal to zero. That is, 

M = Mo, 

A, = Ay - *y (TS + El - E + M p + Pv + p^) - v iPj + mn,», (10) 
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3 Equations and Fluxes 

Here, we derive dynamics equations and thermodynamic fluxes for the system. Neglecting dissipation at this 
point, we assume that the entropy conservation law is valid, 

S + Fu = 0, 

where the entropy flux Fi is determined by 

F = 5(v + w). 

We continue with the equation for the momentum flux found in JtJ , 

ji + U ik:k =0, (11) 

where 

n ik = pw l w k - EiS ik + 2afal (J^p + \f J )) + w ^ k + Wk ^ 1 

( dEi I ( de 
= pWiWk - EiS ik + mn(wiV k + w k Vi) + 2a" af ( — — -s + ( / 



dg a P \ \dg a ' 3 / 
= pwiw k - EiS lk + mn(w t v k + w k Vi) + A ik , (12) 

where we employed the definition (^]) for A^. 

Taking the appropriate equation for the energy flux from , we have 

E + Qi,i = 0, (13) 

where 

/ dH \ up' 
Qi = WiEi + ( e -Q^f) 2~^» + w k U lk . 

To find the second term, we again use the distribution function from (g), 



e (t^- + Vi + Wij f\ = {Vi + Wi ) (ef) + T l{Tx + pv + <j>)^f 
= (Vi + u>i) (ef) - T (pv 91n gp +/) ) = («* + «0 («/> + ^ (ln(l + /)) 
= (Vi + (ef) + T Vl An(l + /) + / In ( ) - Tvi If" pV "' 



/ // V T 

= (uj + Wi) (ef) + TviS - Vi (fe) + v{vP + mn<j) = Wi (ef) + Vi(TS + vP + mn4>). 
For the energy flux, we finally have 

2 

w 

Qi = Wi (ef) + Vi(TS + vP + mncp) —(pWi + mnvi) + w k (pwiW k + mn(wiV k + w k Vi) + A ik ) 

o 

W 

= Wi (ef) + Vi(TS + vP + mruj)) + — {pw t + mnvi) + w k (mnw t v k + A lk ) 



(ef) + wj + p— + Vi ITS + vP + mn U + — + w k A 



= Vi l^TS + vP + mn (0 + ^))+ w i( E - E i) + ( 14 ) 

This formula completes the list of the conventional elasticity theory equations. An additional equation is 
required to govern the quasiparticle degrees of freedom. We now find the time derivative of Pi. We multiply 
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Boltzmann equation (Q) by pi and integrate over the momentum space. We temporarily neglect Umklapp 
processes, which are supposedly rare. If needed, dissipation can be explicitly introduced into the final result. 
In other words, the quasimomentum p is conserved in (normal) collisions, and the term involving the collision 
integral is therefore zero. The left-hand side of the Boltzmann equation gives 



P 



df dH _ df dH\ 
dp dr dr dp J 



-f 



— (— + (■ 
dp \dr 

( de . . dw 



= - / 



dc 



fPi 



d 2 e 



dxi drdp 
The first term can be transformed as 



-fp 

dfck 
' dr dp 



dwA 


_o_i 


f de 


dr J 


dr 1 


K dp~ 


d 2 e 


dw\ 





+ 



\ drdp dr 



dw 

dxi 



df de df 
dr dp dr 



dw dw df 

P"^— / +Pi-£-f +Pi— W 



dxi, 



dr 



dr 



de 
dxi 



fPi 



d 2 



drdp 

= (fp. 



Pi 



df_de_ 
dr dp 



drdp dxi 



. p .<*l( T <b. 

1 dr \ dp 



+ f °L 

<9r<9p dxi 



df 
dr 



df dz ( de 



dz dp \ dr 



dv k 



dT 



df , ,a<t> 

dr dxi 



fp 



dxi 



fz 



dT 

dxi 



P k ~x TO « z ~a~~ 

or dr dr 



dv dzdT 

P^-pT + PiJa — t~ 
dr dp dr 



dcj> _dv d , n . I , dz\ dT 



dS ^dv d , , 

dxi dxi 



dxi 



dr 

dxi 



dp k / dx 

5 lk f\n l -±l-p~^-Mf + l) 

f dp k dz 



| (P(V) + ( /h l±/ +h(/ + 1) )g 

d(j> Jv 
= nm- 1- P-— 

OXi dxi 



dT 

dx k 



dT 



d~r^ + S dx 



Consequently, 



P 



dw 

r 

dxi 



d\ 



d 



dxi dr 



PjW — nm 



dxi 



dxi dr 1 dxi 



nm 



dx-. 



ur 
~~2 



s 



dT 

dxi 



dr 



(i*(w + v)). (15) 



The desirable complete set of the two- velocity elasticity theory equations consists of Eqs. (|l5|), ([l3|) (with Q 
defined by (pi)), and dO) (with Ely defined by (O)). 



4 Boundary Conditions 

We now turn to boundary conditions. They immediately follow from the conservation relations to be satisfied 
at the interface. It is much easier to perform all transformations in the frame of reference of the interface 
itself. All the velocities are therefore taken relative to the boundary. Moreover, we simplify the problem 
by restricting it to the one-dimensional case: all fluxes are supposed to be perpendicular to the flat surface; 
we let the z axis run along this direction. Since no curvature is ascribed to the surface, we ignore capillary 
effects. All calculations done here are valid within the linear approximation. Naturally, boundary conditions 
should depend on the type of the boundary and on the type of the media on the other side of the interface. 
We begin with the situation extensively discussed in literature, the solid-liquid interface [JD. Because the 
possibility of the mass flux through the lattice is taken into account, the results are different, however. 
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Figure 2: Solid-liquid boundary: fluxes in one dimension 
4.1 Solid— Liquid Interface 

The liquid on the other side of the interface (being superfluid) is characterized by the chemical potential p L , 
normal and superfluid densities p^ and pg, normal and superfluid velocities v% and v^, temperature T L , 
pressure p L , and the entropy density S L (see Fig. ||) . 



R + S s (v s + w s ) = S L v%, 

(16) 



w b (E b - E b + A b zz ) + v b (T b S b + m b n b <t> b ) = p L (p>s + Pn<) + S T v n , 



v s m s n s + p s w s = p L s v L s + p%v%, 

A zz~ E t -P ■ 

The superscript S indicates that the appropriate quantities refer to the solid. The first equation is the 
entropy growth condition, where R is the surface dissipative function. The last three equations in ( |l6| ) are 
simply the requirements for the energy, mass, and momentum conservation for the surface, respectively. The 
surface dissipative function must be a positive square form. Using ( fl6|) and (|l0|), it can be expressed as 

RT L = v s (S S (T S - T L ) + m s n s ((t) s - p L )) + w s (E s - E? + A s zz - p s p L - T L S S ) 

= v s (m s n s (<t> 3 - p L ) + S S (T S - T L )) + w s (X s zz + p s {4> S ~ H L ) + S S (T S - T L )) . 

We now recall that the solid-liquid boundary can be either atomically-rough or atomically-smooth, 
depending on the temperature. The nature of the surface may (or may not) impose certain restrictions on 
the dynamics. For both types of the surface, the equation 



A s zz -E?=p L (17) 



is satisfied. 



4.1.1 Rough Surface 

Employing the Onsager principle, we obtain 

v s = a (m s n s {4> s - p L ) + S s \T S - T L )) + r, (\ s zz + p s (0 s - p L ) + S S {T S - T L )) 
w s = 7] (m s n s (^ s - p L ) + S S (T S - T L )) + v (A* + p s (0 s - p L ) + S S (T S - T L )) . 

^ ^ I is positively definite. 

4.1.2 Smooth Surface 

A smooth surface implies immobility of the interface relative to the lattice. That is, 

w s = 0. 



(18) 
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For the quasiparticle gas velocity, we then obtain a restricted version of (18) 

V s = a (m s n s ((t) s - fi L ) + S S (T S - T L )) , 
with the kinetic coefficient a > 0. 



(19) 



4.2 Solid-Wall Boundary 
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Figure 3: Solid-wall boundary: fluxes in one dimension 



By a wall, we imply a macroscopically flat structureless medium, in short, "concrete". The "Solid-wall" 
boundary occurs between solid helium and some normal rigid solid (silica in experiment Q). A "concrete" 
wall is characterized by no mass flux in it (i.e., through the interface). The wall can supply an arbitrary 
energy flux; we let Q denote the flux and T w the wall temperature (see Fig. ||). Concrete is characterized 
by fewer variables than liquid, and the appropriate equations are therefore somewhat simpler. 

Just like for the solid-liquid interface, the actual boundary conditions must depend on the microscopic 
pattern of the surface. One can easily imagine a smooth basal plane of the crystal adjacent to an atomically 
flat concrete wall. This plane must stay at rest with respect to the wall since its motion presupposes the 
creation of new atomic layer nuclei. The plane is similar to the smooth solid-liquid interface, and we can 
naturally say that such an interface is smooth. The boundary condition is then given by w s = v s = 0. 

Another, much more interesting scenario is realized if the interface is slightly tilted with respect to the 
basal plane. Such planes may move by growing additional nodes at the edge. This means that no restrictions 
are imposed on the lattice velocity near the interface. In other words, vacancies are allowed to freely appear 
and vanish on the surface (in this sense, the surface is similar to a grid of dislocations arranged at the 
boundary of the crystal that serve as sources or sinks for vacancies; similar speculations may be found in Q 
in explaining polycrystal plasticity). We call this type of the interface "rough". In this sense, the solid-wall 
boundary can be either smooth or rough. The suggested growth mechanism can be applied only to the rough 
boundary. 



4.2.1 Rough Boundary 

Assuming the boundary to be rough and using the same approach as for the liquid, we write the conservation 
laws 

S s (v s +w s ) =R+Q/T W , 
w s (E s - Ef + Af z ) + v s (T s S s + m s n s (f) s ) = Q, 

v s m s n s + p s w s = 0. 
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Again, R is the surface dissipative function, 

S -fL {E s _ E f + A s z _ p s^s _ T w sS) + s s {T w _ tS) 



RT W = v s 



m n 



S 



A* + 1 - 



9 S 1 

m n 



P \ P 

It must be positive, and for the quasiparticle velocity on the surface we therefore have 



v s = (3 



S S / s s 



xt + i - 



S (T — T ) , 



(20) 



where (3 > is the surface kinetic coefficient. 



5 The Growth Rate 

We now use the equations and boundary conditions obtained above. The physical system discussed in what 
follows is solid helium with elementary excitations represented by phonons and vacancies. We first introduce 
a certain amount of friction between the quasiparticle gas and the lattice. To obtain a physically sound 
result, we again restrict our analysis to one dimension. Furthermore, for simplicity, all our calculations are 
performed within the linear approximation. We can write the quasimomentum density (with the superscript 
S omitted) as K = pkv, where 



PK 



h T 4 



(21) 



Here, Qd is the Debye temperature, c is the velocity of sound in the crystal, and a is the lattice period. 
The last equation is quite obvious. It follows from the fact that in the low-temperature region, the quasi- 
momentum is mainly associated with phonons (the number of vacancies is exponentially small) . The result 
therefore coincides with the one for the mass density (and the momentum density) of the normal component 
of the superfluid, p K - p n ~ T A /h 3 c 5 (see pj). 
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Figure 4: Crystal growth in one dimension 

To describe (rare) Umklapp events, we introduce the appropriate relaxation time parameter tjj. It is a 
between-Umklapp-collision time". From ([l5]), we then have 

= K = -7rmV0 - SWT - — = -mnV0 - SVT - p K — . 
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It is worth mentioning that 777 may well depend on both phonons and vacancies, despite the fact that the 
population of vacancies is far lower than that of phonons. For instance, if the vacancy energy band is 
sufficiently narrow, then the probability of Umklapp processes is significantly higher for vacancies than for 
phonons. This might overcome the low concentration of vacancies. 

Interestingly enough, these formulas allow us to obtain the growth rate for a smooth surface. The 
quasiparticles playing the crucial role here (that of mass carriers) are vacancies, with their mass given by 
m = -mue- 

To estimate the rate, we write the temperature gradient as VT = (T 2 — T{)/h, where the subscripts 1 and 
2 stand for the solid-wall and the solid-liquid interfaces respectively (see Fig.j|). Likewise, for the chemical 
potential we write V</> = (fa — fa)/h. We now use the boundary conditions (119) and (EOh, 



mn\ 



(a + (3)v = af3 \mn(fa - ^) + S(T 2 - T L ) + ^Kzi + I 1 - J S(T W - Ti) 

= a(3^mn(fa - fa) + S(T 2 - T\) + S(T W - T L ) + mn (fa ~ ^ + ~ p ~ S ( tW ~ } 

= ap{-p K h^- + S(T W - T L ) + mn Q (\ zzl - S(T W - Tij) + fa - p L \ }. 

In other words, 

v (± + ~ + P -^\ = S(T W - T L ) + mn fl (\ zzl - S(T W - T,)) +fa- A , 

where (using Eqs. ( |l0| ) and dl7|)) 

A 2Z i = A 2zi - TxSi - En + Ei - fapi = p L - T1S1 +Ei- fap\. 

In equilibrium, \ zz i — 0, fa = p L , and T w = T L = T\. If the liquid temperature and pressure change by 
AT and Ap, we can write an equation for the growth rate v g , 

- Vg ^a + UBEh)=-SAT 
mn \a p tjj J 

+ mn Q (Ap - SATx - TASi + AE ± - pAfa - fa\ Pl + SAT X ) + Afa + S * ^ 

= -SAT + ran (2t + «^E) = - ( S - AT + ,„„ (I - Ap, (22) 

where we used Eq. (^), the thermodynamic equality dp = (dp — SdT)/p for the liquid, and the obvious 
relation v g — —vmnj p. 

We now consider this equation with the second term in the right-hand side equal to zero. This is a usual 
scenario for heat conductivity measurements. The heat flux Q = vTS can then be expressed as 

AT 



Rk i + Rk2 + H/k' 

where Rk i and Rki are the Kapitza thermal resistances on the solid-wall and solid-liquid boundaries and 
k is the heat conductivity of the crystal. Taking the inequalities Rk2 *C Rki, win P, and p — p L <C p into 
account, we immediately obtain 

R 1 TUTS2 



(3TS 27 pk ' 

As a result, the growth rate is given by 

mn ( cat , P-P L 



Strictly speaking, the last equality implies that thermodynamic properties of the crystal mainly depend on 
phonons, while the contribution of vacancies to the effect under consideration is limited to the mass transfer. 
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The growth rate here depends on the overpressure as well as on the temperature difference between the 
liquid and the wall. 

In the real experiment Q], the temperature is lower than 100 mK. In this region, the phonon free path 
is much longer than the experimental cell size and the impact of phonons on the vacancy behavior is 
proportional to a high power of the small ratio T/Qp. Consequently, as Qd oo, we can consider the 
vacancy gas as an independent component and neglect the influence of phonons on it. The crystal growth is 
accounted for by the presence of vacancies; hence, to estimate the growth rate in the experimental situation, 
we can simply substitute the vacancy-only quantities for all variables in Eq. (^2|). Since it seems that there 
were essentially no temperature gradients in the experiment, we consider the isothermal case AT = 0. 

The kinetic coefficients a and [3 arc determined by the vacancy annihilation probabilities in vacancy- 
surface collisions. The vacancy gas velocity v at the interface should be expressed in terms of the accom- 
modation coefficient W (which is the ratio of the number of annihilated vacancies to the total number of 
incident vacancies) as 

v ~ ^j-V T W, (23) 

where Af ~ JmAp/T is the difference between the incident and reflected distribution functions and Vr ~ 
y/T/m* is the thermal velocity. Here, m* is the effective mass near the bottom of the vacancy energy band. 
This mass can be estimated from the energy band width A as m* ~ h 2 /(a 2 A). 

The accommodation coefficient W, like any other inelastic process probability in quasiparticle— surface 
interactions frO] , is approximately the squared ratio of the lattice period to de Broglie wavelength, 

We can now compare Eqs. (p3|), (p^), and (53). For the coefficients, this yields, 

r a I? 

An estimate of the relaxation time tn characterizing the normal (non-Umklapp) vacancy collisions can 
be obtained from 

1 

TN 



naVr ' 



where a ~ a 2 is the vacancy-vacancy scattering cross-section. The Umklapp relaxation time is exponentially 
longer tjj ~ Tpj exp(Aj//T), where Ajj < A is a certain energy specific to the vacancy Umklapp processes. 
Using the obvious relation p K = m*n, we proceed to the growth rate. From Eq. p^), it follows that 



m 2 n 2 (p-p L ) ^ m 2 n 2 (p~p L ) ^ 



(1/q + 1/(3 + m *n/Tu)p 3 (^^/AjTnh/a + hm*n 2 a 2 V T exp(-A c/ /T)) p 3 



-l 



Here, as an estimate, we set pa 3 ~ m and na 3 ~ exp(— eq/T), where eo is the bottom of the vacancy energy 
band. For the facet mobility pf± — v g /Ap introduced in J2|, we have 

a 4 p~ P L ( /A /T h It Ai , 



6 Conclusion 




^- Y—r- Vi eeo/J +~\ a^ u/j • ( 24 ) 



Formula (|2J) provides a reasonable correspondence between the theory proposed here and the experiment 
It suggests three main predictions to be verified in further experiments. 
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• The facet mobility has a maximum at some finite temperature. Should the temperature decrease below 
the point of the maximum, the growth rate will also decrease. Otherwise, if the mobility does not tend 
to zero as the temperature tends to zero, this should be considered as an indication of the presence of 
zero-point vacancies (see 

• The observed growth rate depends on the height of the sample. 

• The crystal grows at the boundary between the solid and the wall. This fact can potentially be observed 
experimentally using some small foreign object frozen into the crystal in its upper part. 
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